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ABSTRACT 

Precision Doppler velocity measurements from the Anglo-Australian Tele- 
scope reveal a planet with a 9.4±0.4 year period orbiting the Ml. 5 dwarf GJ832. 
Within measurement uncertainty the orbit is circular, and the minimum mass 
(msinz) of the planet is 0.64±0.06 Mjup. GJ832 appears to be depleted in met- 
als by at least 50% relative to the Sun, as are a significant fraction of the M 
dwarfs known to host exoplanets. GJ832 adds another Jupiter-mass planet to 
the known census of M dwarf exoplanets, which currently includes a significant 
number of Neptune-mass planets. G J 832 is an excellent candidate for astromet- 
ric orbit determination with asinz = 0.95 mas. GJ832b has the second largest 
angular distance from its star among radial velocity detected exoplanets (0.69 
arc sec) making it a potentially interesting target for future direct detection. 
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1. Introduction 

Most of the known exoplanets orbit late-F, G, or early-K dwarfs, witli masses ranging 
from 0.7 to 1.2 Mq. There are nearly 2,000 such stars within 50 pc brighter than V=8. In 
contrast there are just a handful of M dwarfs brighter than V=8. 

Although M dwarfs make up 70% of nearby stars, their faintness in the optical makes 
them difficult targets for which to obtain precision Doppler velocities. As a result they make 
up as little as 5% of the current planet search targets, and only 11 exoplanets have been 
found to date, orbiting a total of 7 M dwarfs. Just over half of these M dwarf planets are 
less massive than Nept une, leading to spec ulation that M dwarfs typically host either fewer 



planets than G dwarfs (IJohnson et al. 



proto-planetary disks (iLaughlin et al 



20071) , or lower raass p l anets as a result of their sm aller 



2004J : llda fc Linll2005l : iKennedv fc KenvonI 120081 ) 



The primary parameters in planet formation theory are the mass of the central star, 
the mass of the protoplanetary disk, and the metallicity of the system. While it is now 
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), it has been harder to establish the importance of 



stellar mass on planet formation since most of the stars under survey lie in the relatively 
narrow mass range encompassed by late-F, G and K dwarfs. 

The searches which are being done for planets orbiting M dwarfs, will ultimately provide 
the data needed to see if the metallicity-planet relation extends down to the M dwarf regime, 
and whether the mass distribution of exoplanets formed around M dwarfs is similar to, or 
different than, that for more massive host stars. M dwarf Doppler surveys, therefore, have 
the power to address some of the most important questions in exoplanetary science, as they 
extend the mass range of potential exoplanet host stars down to 0.3 Mq. 

We report here a new extrasolar planet in a long period orbit with eccentricity consistent 
with zero, discovered by the Anglo-Australian Planet Search (AAPS). The AAPS program 
is described in Section 2. The characteristics of the host star and our Doppler measurements 
are presented in Section 3. A discussion follows. 



2. The Anglo- Australian Planet Search 



The AAPS began in 1998 January, and is currently surveying 250 stars. Thirty exo- 
planets with msinz ranging from 0.17 to lOMjup have first been discovered by the AAPS 
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Doppler measurements are made with the UCLES echelle spectrometer (jPiego et al.lll990l ) on 
the 3.9m Anglo- Australian Telescope (AAT). An iodine absorption cell provides wavelength 
calibration from 5000 to 6200 A. The spectrometer point-spread function and wavelength 
calib ration are derived f rom the iodine abs option lines embedded on every spectrum by the 
cell JValenti et al.lll995l : ISutler et al.lll996h . Our observing and analysis system has demon- 
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GJ 832 



At 4.93 PC, GJ832 (^8 3 865, HD 204961, HIP 106440) is amongst the nearest stars in 
the sky (jPerryman et al.lll997l ). It is an Ml. 5 dwarf with an optical absolute magnitude and 
colors of My = 10.19, V = 8.66, and B — V = 1.52, and an infrared absolute magnitude 
and colors of Mk = 6.03, K = 4.5 0, and V — K = 4.16. Both Hipparcos and ground- 



based photometry (IKoen et al.l 12002 ) find G J 832 to be photometrically stable at the several 
milli-magnitude level. Gautier et al. ( 2007 ) have combined Nstar^ visible photometry with 
Spitzer far-infrared photometry, to estimate an "infrared flux method" effective temperature 
of 3657 K for GJ832. The Spitzer observations reveal no evidence of mid- or far-infra red 
excess. The radius of GJ832 is estimated to be 0.48 R© (jPasinetti-Fracassini et al.l 120011 ). 



Although accurate metallicities for M dwarfs are problematic, GJ832 is likely to be 
rather metal poor. Matching synthetic spectra to high-resolution spectra of the FeH band 
near 9900A, Schiavon et al. (1997) estimate a metallicity for GJ832 of [Fe/H] = —0.7, and 
a surfa ce gravity of log^r = 4.7. The photometric metallicity calibration of iBonfils et al. 
(l2005al ) gives an estimated metalhcity of [Fe/H] = —0.31 ± 0.2. 



Due to its late spectral type, GJ 832 has not (to date) been subject to detailed spectro- 
scopic analysis, and so to estimate its mass we must rely on either theoretical isochrones, or 
empirical mass-luminosity calibration. The latter indicate a mass for GJ 832 of 0.45±0.05 Mq 
(with the mass unc ertainty being largely due to the scatter about the mass-lumi nosity calibra- 



tion r elationship of iDelfosse et al.l (119981 )). The Padova theoretical isochrones (IMarigo et al. 



20081 ) predict Mk ranging from 5.97 (at 10^ yr) to 6.03 (at 10^° yr) for a 0.45 M© dwarf with 
[Fe/H] = —0.3 which is consistent with the observed luminosity of GJ832. At [Fe/H] = —0.7 
they predict Mk in the range 5.92 (at 10^ yr) to 5.86 (at 10^° yr). Given the difficulty in 
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determining metallicities for M dwarfs, we therefore derive a mass estimate for the primary 
of 0.45±0.05 Mq. 

GJ832 is c hromospherica ll y quie scent. Based on high resolution spectroscopy of the 
Call H&K hues, iTinney et al.l (l2002bl ) report logRl-^K — ~5.10. This would suggest a jitter 
of 3.9 ms~^ using the B — V, My fc Tpff in the most recent stellar "jitter" calibration of 
J.Wright (priv.comm). iBonfils et al.l (l2005bl ) estimate the stellar jitter of GJ832 to be less 
than 2ms^^. GJ832 is among the fainter stars on the AAT program. The signal-to-noise 
of these observations range from 46 to 150 per spectral pixel, with a median of 98, which 
is lower than typical for AAPS targets. Four late dwarfs from the long term AAT program 
are shown in Figure [TJ These stars are shown in order of descending B— V. GJ 887 is an 
especially close match to G J 832 in B-V colour and V magnitude. Based on this we estimate 
the combined velocity uncertainty due to photon statistics, jitter, unknown planets, and 
systematic errors is 5ms~^ for late-K and M dwarfs in the AAPS. This is comparable to 
that estimat e d for late-K and M dwarfs in the Keck program, as shown in Figures 2-4 of 
Butler et all J2008h . 



A total of 32 precision Doppler measurements of GJ832 spanning 9.6 years are listed 
in Tableland shown in Figure [2] (upper panel). The root-mean-square (RMS) scatter of 
the residuals about the mean velocity of this dat a set is 11.6 ms~^. Usi ng the 2-Dimensional 
Keplerian Lomb-Scargle (2DKLS) peridogram of lO'Toole et al.l ( 120071 ) to identify an initial 
period and eccentricity, the subsequent best-fit Keplerian to all 32 epochs of data reduces this 
to an RMS of 5.5 ms~^, and gives a reduced of 1.54 (see Table[2]- a stellar jitter of 3.9 ms~^ 
was used, together with the internal velocity measurement uncertainty for each epoch in 
Table [H to determine reduced x^). These fit parameters strongly suggest the presence of an 
exoplanet with minimum mass msmi of 0.64 M©, period 9.4±0.4yr, eccentricity 0.12±0.11 
(which we consider to be consistent with zero e ccentricity, part i cularl y when the bias against 
measuring zero eccentricities demonstrated by lO'Toole et al.l ( l2008l ) is taken into account) 
and semi-major axis 3.4±0.4AU. 

We have determined the False Alarm Probability (FAP, i.e. the probability that we 
have falsely identified an exoplanet that is not present) for this orbit determ i natio n using 
the Monte Carlo "scrambled velocities" approach described by iMarcy et al.l (120051 ). This 
method tests the hypothesis that no planet is present and the Keplerian fit could have been 
obtained from mere noise, by generating randomly scrambled data sets in which the order of 
velocities are changed but the times remain the same. These are then subjected to the same 
analysis as our actual data set (i.e. identfying the strongest peak in the 2DKLS followed by 
a full Keplerian fit). In this case 2002 random trials were carried out and only one of these 
yielded a less than the value of 1.54 obtained with the original data. The histogram of xi 
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is shown in Figure [3l These results imply a FAP of 0.05% for the GJ 832 planet detection. 



Discussion 



GJ832 at a distance of 4.93pc is one of the nearest known exoplanetary systems. The 
combination of the small distance and relatively long period gives a large angular distance 
from the star of 0.69 arc seconds for an edge-on circular orbit. This is exceeded only by 
e Eri among radial velocity detected exoplanets, and only six other systems exceed 0.2 arc 
sec. G J 832b is therefore a potentially interesting target for direct det ection, although the 
high constrast with the star (likely to be < 10~^; iBurrows et al.l |2004| ) still makes this an 
extremely challenging observation. 

GJ832 is an excellent candidate for astrometric orbit determination. The astrometric 
orbit semimajor axis is asini = 0.95 mas, which is c o mpar able to that of e Eri for which an 



astrometric orbit was determined by lBenedict et al 



J2006h and larger th an that of GJ 876 



which also has an astrometric orbit determination (IBenedict et al.ll2002l ). The astrometric 
orbit would enable the inclination to be determined, removing the current sini uncertainty 
on the mass. 

Seven M dwarfs (including GJ832) are currently known to host as many as 11 exo- 
planets, and these are listed in Table [3] (see table notes for references). As noted earlier, 
determining the metallicities of M dwarfs is notoriously difficult - published metallicity esti- 
mates are available for several of the known exoplanet host M dwarfs, and these are listed in 
the Table. In addition, we have al so derived for all seve n M-dwarfs a photometric metallicity 
estimate, using the technique of iBonfils et al.l (l2005a[). which has the advantage of being 
uniform over all these M dwarfs. On average the ISchiavon et al.l ( 119971 ) metallicities appear 
to b e systematically .3-0.4 dex lower than those derived from the Bonfils et al. calibration. 
The iBean et al.l (120061 ) metallicities are similarly on the metal-poor side of the Bonfils et al. 
results, though not by as much (?» 0.2 dex). In general, the metallicity trends are similar 
across all three calibrations, and it is clear there is a metallicity spread across the observed 
M dwarf exoplanet hosts. 

Based on these metallicity estimates, it would appear that four of the current M dwarf 
exoplanet host stars are somewhat metal-poor, two have about solar metallicity, and one 
is slightly metal rich. Given the well known correlation between stellar metallicity and 
observed exoplanet frequency for F, G, and K dwarf host stars, this metallicity distribution 
for M dwarf host stars is quite unexpected. Whhile the numbers of systems are small there 
is no obvious difference in metallicity between the stars hosting Jupiter-mass planets and 



- 6 - 



those hosting Neptune-mass planets. 



The correlation between high stellar metallicity and planets for late-F, G, and early-K 
dwarfs points toward the core accretion model for planet formation. But there does not ap- 
pear to be strong evidence to date that M dwarf planet formation is strongly correlated with 
high metallicity. This is puzzling, particularly in view of the fact that M dwarfs probably 
have lower mass protoplanetary disks, and therefore would need even higher metallicity than 
a F, G or K dwarf to provide enough solid material (silicates and ices) to build a plane- 
tary core. Obviously it must b e kept in min d that m easuring metallicities for M dwarfs is 
problematic, and that even the iBonfils et al.l (l2005al ) calibration (though empirically based 
and moderately robust) is only good to ±0.2 dex. Nonetheless it is interesting to consider 
possible means by which M dwarf exoplanets could be formed in such a manner as to not 
display the strong metallicity correlation seen in FGK dwarfs. One initially attractive ex- 
planeation is that since M dwarfs are essentially immortal on a Hubble time scale, the vast 
majority of nearby M dwarfs could be old metal poor stars. Unfortunately such an explana- 
tion would a ppear unlikely. The study of M dwa r f kinematics has an extensive and venerable 
history (e.g. I WielenI 119771 : IWeis fc Upgrenlll995l : iReid et al.l 1 19951 1 which has contributed to 
the creation of extensive and sophisticated models of the stellar popul ations present in the 
Solar Neighbourhood (e.g. the Besancon models of iRobin et al.ll2003l ). More recently, the 



availability of huge numbers of M dwarf spectra from the SPSS survey hav e enab led sophis- 
ticated tests of the kinematics of the Besancon models by lBochanski et al.l (120071 ). and have 
substantially born out the Besancon model predictions for M dwarfs. Those models indicate 
that dominant solar neighbourhood M dwarfs will be thin disc members with ages almost 
uniformly spread between 0.1 and 10 Gyr. Thick disc M dwarfs (which would indeed be 
expected to have systematically lower me tallicities) will be present at much lower densities 
(around a factor of one twentieth or less; iRobin et al l boosh . and the probability that they 
would make up four of the seven M dwarf exoplanet hosts would seem to be negligibly small. 

An alternative explanation could be that M dwarf planets might form primarily via 
the disk instability mechanism (see e.g. lBossll2008l . and references therein), rather than via 
core accretion, which would make their formation probability more or less independent of 
metallicity. 

Six of the eleven exoplanets known to orbit M dwarfs have minimum masses less than 
O.lMjup. In contrast, only nine planets with msini< O.lMjup have been found among 
the 216 Doppler vel ocity planets with B-V < 1.2 in the "Catalog of Nearby Exoplanets" 
teutler et allboOGbh . With a minimum mass of 0.64±0.06 Mjup, GJ832b is the fifth jovian 
mass planet found orbiting an M dwarf. The most massive M dwarf planet yet found is 
1.93 Mjup . Since massive planets are by far the easiest ones to find, planets of more than 
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2 Mjup orbiting within 3 AU of M dwarfs must be rare, occuring less than around once per 
300 M dwarfs. 
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Table 1. Velocities for GJ 832 
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Table 1 — Continued 
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Table 2. Orbital Solutions for GJ 832 



Parameter Value 



Orbital period P (days) 3416±131 

Velocity semiamplitude K (ms^^) 14.9±1.3 

Eccentricity e 0.12±0.11 
Periastron date (Julian Date-2451000) 211±353 
a; (degrees) 304±38 

msini (Mjup) 0.64±0.06 
semimajor axis asini (AU) 3.4±0.4 
Nobs 32 
RMS (ms-^) 5.5 
Xl 1-54 
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Table 3. Known M dwarf Exoplanet Hosts 



Host 


Type 


Mk 


V- K 


Phot [Fe/H]^ 


[Fe/H]^ 


[Fe/H]= 


msini(Mjup) 


Refs'^ 


GJ832 


Ml. 5 


6.03 


4.16 


-0.31 


-0.7 




0.64 


1 


GJ876 


M4 


6.64 


5.16 


0.02 


-0.4 


-0.12 


0.019,0.619,1.935 


2,8,9 


GJ849 


M3.5 


5.87 


4.83 


0.16 






0.82 


4 


GJ317 


M3.5 


7.26 


4.97 


-0.23 






1.2 


7 


GJ436 


M2.5 


6.02 


4.61 


-0.02 




-0.32 


0.067 


3 


GJ581 


M2.5 


6.85 


4.72 


-0.26 


-0.1 


-0.33 


0.049,0.016,0.026 


5,10 


GJ674 


M3 


6.57 


4.50 


-0.30 






0.035 


6 



^Photomet r ic [Fe/ H] determined using catalogued V, 2MASS Kg and parallax data, with the 



Bonfils et al.l (j2005al ) relation. 



^Metallicity estimates from ISchiavon et al.l (119971 ) 



^Metallicity estimates from lBean et al.l (120061 ) 



M dwarf exoplanet pro p erties from; 1 - this paper; 2 - iDelfosse et al.l (|l998l): 3 -iButler et al 



(l2004h: 4 - iButler et all (l2006a[): 5 



Johnson et al 



Bonfils et al 



J2005bl): 



6 



Bonfils et al. 



(l2007[l: 



J2007h : 8 - iMarcv et al.l J200lh ; 9 - iRivera et all J2005h : 10 - ludrv et all J2007h 
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Fig. 1. — Four stable late K and M dwarfs from the AAT. 
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HD204961 




-20 k 



1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 

Fig. 2. — Doppler velocities for GJ 832 spanning 9.6 yr. The upper panel shows the measured 
velocities with a best-fit Keplerian over-plotted as a dashed line. The residuals to this fit are 
plotted in the lower panel. The Keplerian orbital parameters obtained listed in see Table Ej 
and strongly suggest the presence of a msin z=0.64 Mjup exoplanet. 
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12 3 4 5 

Scrambled Velocity Best— fit 

Fig. 3. — Assessment of the FAP of the Keplerian model for GJ832. The histogram shows 
the values of xl from 2002 trials with randomly scrambled velocities. Only one of these trials 
had xl lower than the value of 1.54 from the original fit, implying a FAP of 0.05% 



